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Within the frame of macroscopic QED in linear, causal media, we study the radiation force of 
Casimir-Polder type acting on an atom which is positioned near dispersing and absorbing magne- 
todielectric bodies and initially prepared in an arbitrary electronic state. It is shown that minimal 
and multipolar coupling lead to essentially the same lowest-order perturbative result for the force 
acting on an atom in an energy eigenstate. To go beyond perturbation theory, the calculations are 
based on the exact center-of-mass equation of motion. For a nondriven atom in the weak-coupling 
regime, the force as a function of time is a superposition of force components that are related to the 
electronic density-matrix elements at a chosen time. Even the force component associated with the 
ground state is not derivable from a potential in the ususal way, because of the position dependence 
of the atomic polarizability. Further, when the atom is initially prepared in a coherent superposition 
of energy eigenstates, then temporally oscillating force components are observed, which are due to 
the interaction of the atom with both electric and magnetic fields. 

PACS numbers: 12.20.-m, 42.50.Vk, 42.50.Nn, 32.70.Jz 



I. INTRODUCTION 



It is well known that in the presence of macroscopic 
bodies an atom in the ground state (or in an excited en- 
ergy eigenstate) is subject to a nonvanishing force — the 
Casimir-Polder (CP) force — that results from the vac- 
uum fluctuations of the electromagnetic field. CP forces 
play an important role in a variety of processes in physical 
chemistry, atom optics, and cavity QED. Moreover, they 
hold the key to a number of potential applications in nan- 
otechnology such as the construction of atomic- force mi- 
croscopes pj] or reflective atom-optical elements [2] . Over 
the years, substantial efforts have been made to improve 
the understanding of CP forces (for reviews, see Ref. [^). 
Measuring CP forces acting on individual particles is a 
challenging task. Since the early observation of the de- 
flection of thermal atomic beams by conducting surfaces 
Q , measurement techniques and precision have been im- 
proving continuously. More recent experiments have been 
performed with atomic beams traversing between paral- 
lel plates 01 . Other methods include transmission grating 
diffraction of molecular beams , atomic quantum reflec- 
tion 7, 8], evanescent- wave atomic mirror techniques [9|, 
and indirect measurements via spectroscopic means [lOj . 
Proposals have been made on improvements of monitor- 
ing the CP interaction by using atomic interferometry 

The theoretical approaches to the problem of determin- 
ing the CP force can be roughly divided into two cate- 
gories. In the flrst, flrst-principle approach explicit fleld 
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quantization is performed and perturbation theory is ap- 
plied to calculate the body- induced atomic energy shift, 
which is regarded as the potential of the force in lowest- 
order perturbation theory [l^ |l3i y, [l5|, [li [itI, 01, 111 ■ 
The calculations have typically been based on macro- 
scopic QED, by beginning with a normal-mode decom- 
position and including the bodies via the well-known 
conditions of continuity at the surfaces of discontinuity. 
Since in such a (noncausal) approach the frequency de- 
pendence of the bodies' response to the fleld cannot be 
properly taken into account, material dispersion and ab- 
sorption are commonly ignored. As has been shown re- 
cently [l^ , the problem does not occur within the frame 
of a generalized quantization scheme that properly takes 
into account a Kramers-Kronig consistent response of 
the bodies to the field. Clearly, the problem can also 
be circumvented in microscopic QED, where the bod- 
ies are treated on a microscopic level by adopting, e.g., 
harmonic-oscillator models [ij]. In the second, semiphe- 
nomenological approach, the problem is circumvented by 
basing the calculations on linear response theory (LRT) , 
without explicitly quantizing the electromagnetic field 
[2ll[2l[2l|2i,[2i,[2i,[23,[2|. In the ansatz for the force, 
either the field quantities or both the field and the atomic 
entities are expressed in terms of correlation functions, 
which in turn are related, via the fluctuation-dissipation 
theorem, to response functions. 

At flrst glance one would expect the result obtained 
from exploiting LRT to be more generally valid than the 
QED result obtained in lowest-order perturbation the- 
ory. In fact, this is not the case. In both approaches, 
it is not the exact atomic polarizability that enters the 
expression for the (ground-state) CP force but the ap- 
proximate expression which is obtained in lowest-order 
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perturbation theory and which effectively corresponds to 
the atomic polarizability in free space. Since the structure 
of the electromagnetic field is changed in the presence of 
macroscopic bodies, the atomic polarizability is expected 
to change as well. It is well known that the atomic level 
shifts and broadenings sensitively depend on the mate- 
rial surroundings. In particular, when an atom is situated 
very close to a body, the effect can be quite significant 
(see, e.g., Refs. [2^, [30]), thereby changing the atomic 
polarizability. As a result, a position-dependent polariz- 
ability is expected to occur, which prevents the CP force 
from being derivable from a potential in the usual way. 



A way to derive a more rigorous expression for the CP 
force is to base the calculations on the exact quantum- 
mechanical center-of-mass equation of motion of the 
atom as we shall do in this paper. The calculations are 
performed for both minimal and multipolar coupling, and 
contact is made with earlier studies of the center-of-mass 
motion of an atom in free space, with special emphasis 
on the so-called Rontgcn interaction term that appears 
in the multipolar Hamiltonian [31i, j52,, 33, 34, 35j. Af- 
ter taking the expectation value with respect to the in- 
ternal (electronic) quantum state of the atom and the 
quantum state of the medium-assisted electromagnetic 
field, the resulting force formula can be used to calcu- 
late the time-dependent force acting on a nondriven or 
driven atom that is initially prepared in an arbitrary (in- 
ternal) quantum state. In this paper, the force formula 
is further evaluated for the case of a nondriven, initially 
arbitrarily prepared atom, by assuming weak atom-field 
coupling treated in Markovian approximation. It is worth 
noting that the theory, being based on the quantized ver- 
sion of the macroscopic Maxwell field, with the bodies 
being described in terms of spatially varying, Kramers- 
Kronig consistent complex permittivities and permeabil- 
ities [3^ |33], also applies to left-handed materials jsl] 
where standard quantization runs into difhculties. 



The paper is organized as follows. After a brief sketch 
of the quantization scheme (Sec. |TT|, in Sec. |TTT] atten- 
tion is focused on the perturbative treatment of the CP 
force acting on an atom in an energy eigenstate, and pre- 
vious results [20] obtained for dielectric surroundings of 
the atom are extended to magnetodielectric surround- 
ings, including left-handed materials. In Sec. lIVI the exact 
center-of-mass Heisenberg equation of motion of an atom 
and the Lorentz force therein are studied, and Sec. |V] is 
devoted to the calculation of the average force, with spe- 
cial emphasis on a nondriven atom in the weak-coupling 
regime. Finally, a summary and some concluding remarks 
are given in Sec. IVII 



II. SKETCH OF THE QUANTIZATION 
SCHEME 



A. Minimal coupling 

In Coulomb gauge, the minimal-coupling Hamiltonian 
of an atomic system (e.g., an atom or a molecule) consist- 
ing of nonrelativistic charged particles interacting with 
the electromagnetic field in the presence of macroscopic 
dispersing and absorbing bodies reads [l^j 



H = 



E 

A— e,m 



r I duj hujil{r,uj)ix{r,u!) 
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-fi / dVp^(r)^A(r)+ / d-'rpAirMr), (1) 



where 



and 



(^A(r) = / dV 



M(r') 



47reo r — r' 



(2) 



(3) 



are the charge density and scalar potential of the parti- 
cles, respectively. The particle labeled a has charge 
mass TOq, position r^, and canonically conjugated mo- 
mentum Pq. The fundamental Bosonic fields fA(r, tu) [and 
fj^(r,a;)] which can be related to noise polarization (for A 
— e) and noise magnetization (for X=m), respectively, are 
the dynamical variables for describing the system com- 
posed of the electromagnetic field and the medium in- 
cluding the dissipative system responsible for absorption. 



/A,(r,c.),/|,,,(r',c.')l ^SxySu'd{r-r')d{LO-uj'), (4) 



0. 



(5) 



Note that the first term on the right-hand side of Eq. ^ 
is the energy of that system. Further A(r) and 0{r) are 
the vector and scalar potentials of the medium-assisted 
electromagnetic field, respectively, which in Coulomb 
gauge are expressed in terms of the fundamental fields 
fA(r, w) [and f|(r, w)] as 

A(r) = / da)(iw)"^E-L(r,cj) +H.C., (6) 
Jo 



-V^(r) 



dwE"(r,w) +H.C. 



(7) 



where 



E(r,c.)= ^ 

A— e,m 



dVGA(r,r',c^)fA(r',c^), (8) 
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Ge(r, r', u;)^i^J Ini£(r', u;) G{r, r', u;), (9) 



TTEo 



G„(r,r',Lj) = —lmK{r',uj) 

C V TTEo 



x[G(r,r',c^)xVr-], 



(10) 



with [G(r,r',cj)x Vr'] =ejfe;(9;Gifc(r,r',Lj) and n{v,uj) 

= /i~^(r,a;). Here and in the following, transverse and 
longitudinal vector fields are denoted by _L and ||, re- 
spectively, e.g., 

E-L(ll)(r,cj) = y" dV<5-L(ll)(r-r')E(r',tj), (11) 



with 



and 



4(r) 



1 

V 47rr 



(12) 



(13) 



being the longitudinal and transverse dyadic 5 functions, 
respectively. 

In Eqs. ([9]) and (fTO|. G(r, r', cj) is the (classical) Green 
tensor, which in the case of magnetodielectric matter 
obeys the equation 



V X K(r, cj)V X 



■e(r,cj) 



G(r, r',w) = 5{v-y') 
(14) 



together with the boundary condition 
G(r,r',tj)^0 for|r-r'| 



(15) 



Note that the (relative) permittivity e{v,uj) and perme- 
ability pl{v^lo) of the (inhomogcneous) medium are com- 
plex functions of frequency, whose real and imaginary 
parts satisfy the Kramers-Kronig relations. Since for ab- 
sorbing media we have Ime(r,a') > and Im/z(r, w) > 
^ Im K(r, uj) < 0, the expressions under the square roots 
in Eqs. ^ and pUj) are positive. It should be pointed out 
that the whole space is assumed to be filled with some 
(absorbing) media, in which case the aforementioned con- 
ditions for Im£(r,Ci;) and Im/z(r,w) ensure that the dif- 
ferential equation (fT4|) together with the boundary con- 
dition psp presents a well-defined problem. However, as 
this assumption allows for both e{Y^Lo) and fi{r,uj) to be 
arbitrarily close to unity (i.e., for arbitrarily dilute mat- 
ter), it is naturally possible to include vacuum regions in 
the theory, by performing the limit e(r, lj) — > 1, ^{v,uj) 
^ 1 in these regions after having calculated the desired 
expectation values of the relevant quantities as functions 
of e(r,u;) and ^i{y,uo). 

The Green tensor has the following useful properties 



G(r,r',tj) = G~^{v',v,uj) 



(17) 



j dh[luvK{s,uj)[G{v,s,uj) X Vs] [VsX G*(s,r',c^)] 
+ -^Ime(s,w)G(r,s,w)G*(s,r',w)| = ImG(r,r',w). 



(18) 



Combining Eq. HH) with Eqs. jS]) and ^ yields 
/ d^sGAifc(r,s,w)G^ji,(r',s,a;) 



uP'lm Gij (r, r' ,uj). 



(19) 



Note that in Eq. (fT9|) and throughout the remaining part 
of this paper, summation over repeated vector indices is 
understood. 

The total electric field is given by 



E{r) = E(r) - V^A{r), 



where 



E(r) 



da;E(r,tj) + H.c. 



(20) 



(21) 



with E(r,ti;) from Eq. ([8]). Accordingly, the total induc- 
tion field reads 



where 



i?(r)=B(r)=/ dwB(r,w) -f-H.c, (22) 



(23) 



B(r,w) = (itj)~^V X E(r,w). 



Finally, the displacement and magnetic fields are given 
by 

^(r)=D(r)-eoV^A(r) 



D(r,cj)+H.c. -eoV(^A(r), (24) 

(25) 



H{r)=U{r)= dwH(r,tj) + H.c., 

^0 



where 



D(r,u;) = eoe(r,w)E(r,w) 



— Ime(r, w) fe(r, w). 



(26) 



H(r,[jj) = KoK(r,uj)'B{r,uj) 



ImK(r,w)f„(r,L^). (27) 



Assuming that the atomic system is sufficiently local- 
ized, and introducing shifted particle coordinates 



G*{r,r',u;) = G{r,r\-Lj*), 



(16) 



(28) 
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relative to the center of mass 



(29) 



(niA = X^a'^a)' apply the long- wavelength ap- 

proximation by expanding the fields A(r) and ip{r) 
around the center of mass and keeping only the lead- 
ing nonvanishing terms of the respective field operators. 
For a neutral atomic system, 



9A 



(30) 



this is just the familiar electric dipole approximation, and 
the Hamiltonian ^ simplifies to 



H = Hf + Ha + Haf, 



(31) 



where 



A— e,m 



Hf= / 'i'^ / dc^^fl(r,^)fA(r,c^), (32) 



^A^E|f + i/d''^^A(r)^A(r) 
Haf = dV(^(r)|r=fA 



(33) 



V — PaA(f a) 



(34) 



with 



d ^^qaVa = Yqara (35) 

a a 

being the total electric dipole moment. 

B. Multipolar coupling 

Let us turn to the multipolar coupling scheme widely 
used for studying the interaction of electromagnetic fields 
with atoms and molecules. Just as in standard QED, so 
in the present formalism [36,. ,37.], the multipolar Hamilto- 
nian can be obtained from the minimal-coupling Hamil- 
tonian by means of a Power-Zienau transformation, 



U — exp 



- / d3rPA(r)A(r) 



where 



PA(r) = E^a^a / dA(5(r - fA - Af„). 



(36) 



(37) 



For a neutral atomic system, the multipolar Hamiltonian 
[which is obtained by expressing the Hamiltonian ([1]) in 



terms of the transformed variables] can be given in the 
form of (see Appendix \Q 

A— e.m ^ 



^ / dVPi(r)- / dVPA(r)E'(r) 



E- 



2ma 
where 
B,a{r) = qa&a{r) ^ 

and 

0a(r) =r 



dVH„(r) X B'(r) 



(38) 



m 
niA 



-Y^q^&,ir) + ^pA{r) (39) 



(40) 



dA AJ(r — f A — Ar^). 



Note that due to the unitarity of the transformation (j36|) . 
the transformed variables of the atomic system f ^ — f „ 
and and the transformed field variables i')^{r,uj) and 
f^^(r, Lu) obey the same commutation relations as the un- 
transformed ones. Needless to say that the transformed 
fields E'(r) and B'(r) are related to the transformed fields 
fj^(r, w) and f^^(r, w) according to Eq. ([HD and Eqs. (^1]) - 
(1^ . with primed quantities instead of the unprimed 
ones. The Hamiltonian ((38|) can be regarded as the gener- 
alization of the multipolar Hamiltonian obtained earlier 
for moving atoms in vacuum [sil, [H, [H, [H, [11] to the 
case where dispersing and absorbing magnetodielectric 
bodies are present. In particular, it can be used to de- 
scribe effects specifically due to the translational motion 
of an atomic system such as Doppler and recoil effects. 

Applying the long-wavelength approximation to the 
fields E'(r) and B'(r) in Eq. ((38)) . which is equivalent to 
approximating S{r — f a — Afq,) by S{r — f a) in Eqs. ([57]) 
and ([iU)) . respectively, i.e., 



PA(r) =d,5(r-fA), 



0Q(r) = irQ(5(r-fA), 



(41) 



(42) 



thus 



Hc.(r) = yja5{r - va) + 7r^dS{r - va), (43) 

2mA 

we obtain the multipolar Hamiltonian in long-wavelength 
approximation. 



H — Hp + H'p^ + H'j^p , 



(44) 
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with 



III. VAN DER WAALS POTENTIAL 



i?^= ^ j d?r r du;nwi'^{v,u:)ii{v,uj), (45) 

A— e,m 

Hip E -dE'(f a) + ^ Pifo X B'(fA) 



mA 



p:^dxB'(fA), 



where 



Pa 



(47) 



(48) 



According to Casimir's and Polder's pioneering con- 
cept , the CP force on an atomic system near macro- 
scopic bodies is commonly regarded as being a conserva- 
tive force. In particular, it is assumed that for an atom in 
an eigenstate \l) of the atomic Hamiltonian the position- 
dependent shift of the corresponding eigenvalue due to 
the (electric-dipole) interaction of the atomic system with 
the body-assisted electromagnetic field is the potential, 
also referred to as van der Waals (vdW) potential, from 
which the CP force can be derived, where the calculations 
are usually performed within the frame of lowest-order 
perturbation theory. In this picture, the center-of-mass 
coordinate is a parameter rather than a dynamical vari- 
able (fa ^ Fa). Following this line, we first extend pre- 
vious results [20] , and show that minimal and multipolar 
coupling schemes yield essentially the same expression for 
the force. 



is the (canonical) momentum of the center of mass, and 



Pa = Pa Pa 

mA 



(49) 



denote shifted momenta of the particles relative to the 
center of mass. The first two terms on the right-hand side 
of Eq. (|47p represent electric and magnetic dipole interac- 
tions, respectively, the next two terms describe the (gen- 
eralized) diamagnetic interaction of the charged particles 
with the medium-assisted electromagnetic fields, while 
the last term describes the Rontgen interaction due to the 
translational motion of the center of mass. In particular, 
in (generalized) electric dipole approximation, Eq. (|47[) 
reads 

H;F = -dE'(fA) + — pkdxB'(rA). (50) 
mA 

Recall that the transformed medium-assisted electric 
field E'(f) is related to the physical one, E(f), accord- 
ing to Eq. jXI)) . 

If the center-of-mass coordinate is treated as a (classi- 
cal) parameter (fAi-+rA), then Eq. ([39]) reduces to 



3Q(r) = ga0a(r). 



(51) 



which corresponds to the limit ma/mA — > 0. Hence 
Eq. (liT)) becomes 



iJ^F = -dE'(rA) + E ^ " X B'(rA) 



(52) 



If the paramagnetic and diamagnetic terms are omitted, 
the interaction Hamiltonian simply reduces to the first 
term on the right-hand side of Eq. (|52p . 



A. Minimal coupling 

We start from the minimal-coupling Hamiltonian in 
electric dipole approximation as given by Eqs. (|3ip - ([Ml) 
together with Eq. ([55]) (fa ^ t^a)- Let \n) denote the 
eigenstates of the multilevel atomic system and write Hx 
[Eq. (IMI)] as 



HA = Y.E„\n){n\. 



(53) 



To calculate the leading-order correction to the unper- 
turbed eigenvalue of a state |/)|{0}) due to the pertur- 
bation Hamiltonian ((34| [|{0}), ground state of the fun- 
damental fields fA (r, uj)] , we first note that the first two 
terms have no diagonal elements. Thus they start to con- 
tribute in second order. 



Wfe; + W 



dV 



k A— e,m 

|(;|({0}|dV^(r)|r=r. -E— P"A(rA) 



x|{l,(r,c.)})|fc) 



(54) 



{V ^ principal part), whereas the third term starts to con- 
tribute in first order, 



Aii?,^(/|({0}|E^A2(rA)|{0}>|0 



2ma 



(55) 



Here, |{lA(r, cj)}) =f|(r, w)|{0}) denotes single-quantum 
Fock states of the fundamental fields, and 

u;u = {Eu-Ei)/h (56) 

are the atomic transition frequencies. Since AiE'; and 
l^2Ei are quadratic in the coupling constant [Eqs. (|B9p 
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and (jB10[) in Appendix[B] , thus being of the same order of 
magnitude, the leading-order correction to the eigenvalue 
is given by 



AEi = AiEi + A2E1 



(57) 



A straightforward but somewhat lengthy calculation 
yields (see Appendix [C)) 



AjjJ 



u;ki + w 
[lmG(rA,rA,w) - Im IIg" (fa, ta, t^)] 
^2lnillGll(rA,rA,c^)}dfei, (58) 



with 



i/fe 



m\k) 



(59) 



being the dipole matrix elements. 

Since the atomic system should be located in a free- 
space region, the Green tensor in this region is a linear 
superposition of the (translationally invariant) vacuum 
Green tensor G^^^ and the scattering Green tensor G^^-* 
that accounts for the spatial variation of the permittivity 
and permeability. 



G(r, r', w) = G(°) (r, r', c^) + G^^^ (r, r',u). 



(60) 



As a consequence, the eigenvalue correction AEi can be 
decomposed into two parts. 



AEi ^ AE. 



(0) 
7 



ae\^\vp:). 



(61) 



The TA-independent term AEf'^ associated with the 
vacuum Green tensor gives rise to the vacuum Lamb 
shift and is not of interest here. The rA-dependent term 
AE^'^\v^)^ associated with the scattering Green tensor, 
is just the vdW potential sought, 

C/,(rA) - Ae\^\vp:) = A^e\^\v^) + A2E\^\v^). (62) 

Hence from Eq. ([55]) [G(rA, ta, t^) >^ G^^^ (rA, ta, w)] we 
derive, on recalling Eq. (jl6p and changing the integration 
variable from —uj to cj. 



Mo 



iOkl+i^ 

[ujki^ [G(i) (rA, TA, Lo) - IIg(i)II (rA, ta, uj)] 
-c.2|lG«ll(rA,rA,c.)}-7'jJ|" 



UJkl - UJ 



{ukiuj [G(i) (rA, TA, Lo) ~ IIg(i)II (rA, ta, lo)] 



11(5(1)11 (- 



rA,rA,t^ 



)} 



(63) 



This equation can be greatly simplified by using con- 
tour-integral techniques. G^^''(rA, ta, i^) is an analytic 



function in the upper half of the complex frequency plane, 
including the real axis (apart from a possible pole at 
Lu ^0). Furthermore, knowing the asymptotic behaviour 
of the Green tensor in the limit oj ^ (cf . Ref. fsT] ) , one 
can verify that all integrands in Eq. (j63p remain finite 
in this limit. We may therefore apply Cauchy's theorem, 
and replace the principal value integral over the posi- 
tive (negative) real half axis by a contour integral along 
the positive imaginary half axis (introducing the purely 
imaginary coordinate to — iu) and along a quarter circle 
with infinitely large radius in the first (second) quad- 
rant of the complex frequency plane plus, in the case of 
LiJik > 0, a contour integral along an infinitesimally small 
half circle around Lu — LUik {lo — —ujik) in the first (second) 
quadrant of the complex frequency plane. The integrals 
along the infinitely large quarter circles vanish due to the 
asymptotic property 



lim ^G(^^(r,r,w) =0 

(cf. Ref. [13), so we finally arrive at 

f/^rA) = ;7r(rA) + t/^i-A), 

where 



Ui (rA) = — 2^/ du 



(64) 



(65) 



X 



dikG'^^\rA,rA,iu)dki (66) 
is the off-rcsonant part of the vdW potential, and 

UlirA) = -^ioJ2QiuJlk)uJfk 

k 

X difeReG(^)(rA,rA,W(/c)dfci (67) 

[8(z), unit step function] is the resonant part arising from 
the contribution from the residua at the poles. Note that 
Ui{r\) vanishes when the atomic system is in the ground 
state. For an atomic system in an excited state, U}'{rA) 
may dominate [/"(fa). 

The CP force can be derived from Eq. (|55|) according 

to 



F^rA) = -VAUiivA) 



(68) 



CVa = ^a)- a formula of the type of Eq. (ISsp t ogether 
with Eqs. ([TO]) and ([^7)1 was first given in Ref. [23| within 
the frame of LRT. 

To give Eq. (j66p in a more compact form, we introduce 
the generalized atomic polarizability tensor 



amn(t^) 



-y 



U!kn - LO - i{rk + Tm)/2 
dkn <X' dfnk 

LOkm +1^ + i{Tk + r„)/2 



(69) 
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where Hikm are the shifted (renormahzed) transition fre- 
quencies and T], are the excited-state widths. Following 
Ref. [3^, we may regard 

oli{lo) ^ a.ii{Lo) (70) 

as being the ordinary (Kramers-Kronig-consistent) po- 
larizability tensor of an atom in state |/). Hence we may 
rewrite Eq. (|66|) as 

(71) 

where 

aj°) {lo) = lim ^ V , . dik ® dki (72) 

is the polarizability tensor in lowest-order perturbation 
theory, which can be obtained from Eq. ([70]) together 
with Eq. (f69|) by ignoring both the level shifts and broad- 
enings. In particular for an atom in a spherically symmet- 
ric state, we have 

a(°)(.) = = lim A ^ , . Ml 

k kl 

(73) 

{I, unit tensor), so that Eq. (|7T|) reduces to 

C/r(rA) - ^ y'"diiu2a|°)(zii)TrG(i)(rA,rA,»u), 

(74) 

and Eq. ([67|l simplifies to 

;7/-(rA) - -^5]e(c.,,)c.f,|d,,p 

xTr[ReG(i)(rA,rA,Wifc)]. (75) 

Note that 

ai{iu) ~ OLf\iu) (76) 

is typically valid for an atomic system in free space, be- 
cause of the smallness of the level shifts and broadenings 
that result from the interaction of the atomic system with 
the vacuum electromagnetic field. 

Equation ([55]) together with Eqs. (pS)) and ([57]) can be 
regarded as being the natural extension of the QED re- 
sults obtained on the basis of the familiar normal-mode 
formalism, which ignores material absorption. Moreover, 
it does not only apply to arbitrary causal dielectric bod- 
ies, but, to our knowledge, it first proves applicable to 
magnetodielectric matter such as left-handed material, 
for which standard quantization concepts run into diffi- 
culties. Note that all information about the electric and 
magnetic properties of the matter is contained in the 
scattering Green tensor. 



Finally, let us briefly comment on the ground-state po- 
tential as given by Eq. ([7T|) for / = 0. In terms of an inte- 
gral along the positive frequency axis, it reads 

C^o(rA) = — TT- I dww 
Jo 

X Im{Tr[4")(c.)G«(rA,rA,c.)]}. (77) 

An expression of this type can also be obtained by us- 
ing the methods of LRT [H, [2^ . It allows for a simple 
physical interpretation for the ground-state CP force as 
being due to correlations of the fluctuating electromag- 
netic field with the corresponding induced electric dipole 
of the atomic system plus the correlations of the fluctuat- 
in^ electric dipole moment with its induced electric field 



B. Multipolar coupling 

Let us now consider the multipolar Hamiltonian in 
long-wavelength approximation as given by Eqs. (|44p - 
gSl) together with Eq. and write H'^ [Eq. dH])] in 

the form of Eq. (|53p . In contrast to the electric dipole ap- 
proximation considered in the minimal coupling scheme, 
the present Hamiltonian also includes magnetic interac- 
tions. One might therefore expect that the leading-order 
corrections to the unperturbed eigenvalues are given by 
the second-order corrections due to the dipole interac- 
tions (linear in the field variables) plus the first-order 
correction due to the diamagnetic interaction (quadratic 
in the field variables), all of these contributions being 
quadratic in the coupling constant. However, one can 
show [Eqs. (|BT6| - (|BT5| in Appendix [B] that the second- 
order eigenvalue correction due to magnetic dipole inter- 
action is smaller than that due to the electric dipole inter- 
action by a factor of {ZeSOioY, where ZcS is the effective 
nucleus charge felt by the electrons giving the main con- 
tribution to the energy shift, and ao is the fine-structure 
constant. The current formalism based on Hamiltonian 
(U]) only treats nonrelativistic atomic systems, which are 
characterized by Z^sao <C 1 ^3 ; we can safely neglect 
the correction arising from the magnetic dipole interac- 
tion. Furthermore, the first-order correction arising from 
the diamagnetic term can be shown to be smaller than 
the second-order correction due to the electric dipole in- 
teraction by the same factor (ZoffOfo)^, so we can disre- 
gard it for the same reason. 

In summary, the main contribution to the eigenvalue 
shift of a state |0I{0'}) [|{0'}), ground state of the 
transformed fundamental fields f^(r,w)] is the second- 
order correction due to the electric dipole interaction in 
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Eq. dill), i.e., 

AEi = A2E1 



k X= 



du 



+ w 



X /dV|(Z|({0'}|-dE'(rA)|{l',(r,c.)})|fc)| 



(78) 



[\{l'^{r,uj)}) = {'^ {r,uj)\{0'})]. After some algebra it can 
be found that (see Appendix [Cjl 



duj 

^fci 



- UJ 



X d;fcImG(rA,rA,t^)dfe; 



(79) 



We now apply the same procedure as in Sec. lIII Al below 
Eq. ([55)1 . Replacing the Green tensor by its scattering 
part and transforming the frequency integral to imag- 
inary frequencies using contour integral techniques, we 
arrive at exactly the same form of the vdW potential as 
given in Eq. together with Eqs. ([M)) and ((57)) . It 
is worth noting that the two schemes lead to equivalent 
results only if in the minimal-coupling scheme the 
coupling term is properly taken into account. 



IV. CENTER-OF-MASS MOTION AND 
LORENTZ FORCE 

Atomic quantities that are related to the atom-field in- 
teraction can drastically change when the atomic system 
comes close to a macroscopic body, the spontaneous de- 
cay thus becoming purely radiationless, with decay rates 
and level shifts being inversely proportional to the atom- 
surface separation to the third power [295 . Clearly, in this 
case approximations of the type (|76p cannot be made in 
general and the perturbative approach to the calculation 
of the CP force becomes questionable. Moreover, when 
the atomic system is not in the ground state, then dy- 
namical effects can no longer be disregarded. To go be- 
yond perturbation theory, let us first consider the center- 
of-mass Newtonian equation of motion and the Lorentz 
force therein. 



A. Minimal coupling 

As has been shown 37], the Heisenberg equations of 
motion governed by the minimal-coupling Hamiltonian 



(80) 



lead to the well-known Newtonian equations of motion 
for the individual charged particles, 

TOafa = -B(ra) + 5 raXB{ra)-B{ra)xra >■ (81) 



Summing Eq. ((8T|l over a, recalling definition ([29)) . and 
using Eqs. ([^0]) and together with the relationship 



^9aVa(^A(fa) = 

a 

CVa = ^„), we derive 

mATA = F, 
where the Lorentz force takes the form 



J dV[pA(r)E(r)-|-jA(r) xB(r)], 



(82) 



(83) 



(84) 



with charge density pA (r) and current density ja (r) being 
defined by Eq. © and 



jA(r) ^ ^^Qa [fQ<5(r - Ta) + 5{r - ?„)£■„ 



(85) 



respectively. It can be shown [3l|, \3Q, |4]| that for neu- 
tral atoms the atomic charge and current densities can 
be expressed in terms of atomic polarization and magne- 
tization according to 



PA{r) = -VPA(r) 



(86) 



and 



jA(r) = PA(r) -t- V X MA(r) -f V x MR(r), (87) 
respectively, where 

a 

Mn{r) = i [PA(r) x f-A - i-A x PA(r 



(89) 



with PA(r) and ©^(r) from Eqs. ([371) and (gO]), respec- 
tively. Note that the last term in Eg. (l87|) represents 
the so-called Rontgen current [4l|, |42| |. which is a fea- 
ture of the overall translational motion of any aggregate 
of charges. 

Inspection of Eqs. (gOl), and ^ shows that 
the relations 

V®PA(r) = -VA(8)PA(r), (90) 
V (g> Ma(r) (r) = - Va ® Ma(r) (r) (91) 

{'Va = '^ta) ^re valid. We therefore may write, on recalling 
Maxwell's equations. 



dV [VPA(r)]E(r) 

Va J dV [PA(r)E(r)] -I- J dVPA(r) x ^(r), (92) 
d3r[VxMA(R)(r)] x B(r) 

= VAydV[MA(R)(r)B(r)]. (93) 
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Substituting Eqs. ([86]) and ([STD into Eq. (jM]) and us- 
ing Eqs. (inH) and ([M]) . we may equivalently express the 
Lorentz force as 

F = Va( / dVPA(r)E(r) 



dt 



dV[MA(r)+MR(r)]B(r) 



dVPA(r) X B(r). 



(94) 



In long-wavelength approximation, Eqs. (j88p and 
simplify to [recall Eqs. (jtl]) and ([i^ ] 

MA(r) = i ^ [(^(r ~ fA)rQ x 



-r„ X ra(5(r- ta)] (95) 



and 



MR(r) - i [5 (r - f A)d x f a - x d 5(r - f a)] , (96) 



respectively, so that the Lorentz force ([941) can be written 
as 

F = Va| dE(fA) + 5 X! 9ai'"B(fA) x F„ 
-f- if AB(f a) X d| + [d X B(f a)] . 



(97) 



Further, we calculate 



d_ 

di 



dx B(fA) 



i/,dxB(fA) 



= d X B(fA) + d X B(r) 



-d X i 



r=rA 

FA Va (S) B(f a) + B(f a) (E) Va^a 



(98) 



Comparing the different terms in Eq. (j97|l . one can show 
[Eqs. ((BT9| . (|B20| . and (|B22)) - (|B24l) in Appendix[B] that 
the second term in curly brackets is typically smaller than 
the first one by a factor oi v/c + Zcscto{v, velocity of the 
center of mass), while the third term is smaller than the 
first one by a factor oiv/c. Similarly, we find [Eqs. (|B25[) - 
(|B27)) in Appendix E] that the third term in Eq. ^ 
is smaller than the first two terms by a factor of v/c. 
Thus in the nonrelativistic limit considered throughout 
the current work [cf. Hamiltonian ([T|)] we can set 



V[dE(r)]+-[dxB(r)] 



(99) 



In the absense of magnetodielectric bodies, Eq. 
reduces to earlier results derived within the multipolar 
coupling scheme for an atom interacting with the electro- 
magnetic field in free space [H, [S^. However, it should 
be pointed out that here the electric and magnetic fields 
E(r) and B(r), respectively, are the medium-assisted 



fields as defined by Eqs. ([2T|) and ([22]) [together with 
Eqs. dS]) and Thus Eq. jM]) or, in electric dipole 

approximation, Eq. (|99|1 determine the force acting on 
an atomic system in the very general case of dispersing 
and absorbing magnetodielectric bodies being present — 
a result that has not yet been derived elsewhere. 



B. Multipolar coupling 

Using the multipolar Hamiltonian psp. we obtain, on 
recalling that f ^ = f „ , 



Too 



f-„ = -[i?,TO„f„] ^p'^+ /d3rH„(r)xB'(r). (100) 



Summing Eq. (jlOOp over a and taking into account 
Eqs. (IMl) and yields 



(101) 



H,p'^ + J dVPA(r) X B'(r) 

To evaluate the different contributions to the first term 
in Eq. (fT02|) . we first recall Eq. and note that 



TOArA = Pa + y dVPA(r) x B'(r) 
Equation (jlOip leads to 
TOAf A = F = 



^/d3rVPl(r)=0. 



Further, we derive, on recalling Eq. (jlOOp . 



(103) 



2to 
= -Va 



/" dV i ^ [j„ X H„(r) - H„(r) x f„] B'(r). 

(104) 



Substituting Eqs. (fTOa]) and (fTOi]) into Eq. (fTOS]) . with H 
as given in Eq. ()38|) . we eventually obtain 



F = Va^ / dVPA(r)E'(r) 



^y"dV ^ [Ha(r) X f„ -f-„ X H„(r)] B'(r) 

(105) 



^ / dVPA(r)xB'(r) 



It can be shown (see Appendix [D| that Eq. (|105p is iden- 
tical to Eq. (|M)) . 

It is not diflicult to see [recall Eqs. (|3T|) and ([33])] that 
in long-wavelength approximation Eq. p05p takes the 
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form of Eq. ^7^, but with E'(fA) and B'(fA) in place 
of E(fA) and B(rA), respectively. The time derivative 
d[d X B'(fA)]/dt can then be calculated to give an ex- 
pression of the form of Eq. (^5)) with B(fA) replaced 
by B'(fA)- Obviously, in the nonrelativistic limit we are 
left with an expression similar to Eq. ((99|) . It should 
be pointed out that Eqs. (HH) and (HH) with E(fA) 
and B(fA) replaced by E'(fA) and B'(rA), respectively, 
yield exactly the same force as the equations with the 
unprimed quantities, although the physical meaning of 
E'(fA) is different from that of E(fA) [recall that B'(fA) 

= B(fA)]. 

It is worth noting that the results of this section can 
serve as an example to illustrate that the electric dipole 
approximation has to be employed with great care. If 
in electric dipole approximation the Rontgen interac- 
tion primarily related to the induction field had been 
disregarded and Eq. (|50|1 without the second term on 
the right-hand side had been used, then in the result- 
ing expression for the force the time-derivative term, i.e., 
the magnetic part of the force, would have been lost. 
Note that the pressure exerted by external laser fields on 
macroscopic bodies can be dominated by this magnetic 
force [H, HJ, which contrasts with arguments [33 . |45| 
that the contribution of this term to the radiation force 
on atoms can be neglected. 

V. AVERAGE LORENTZ FORCE 

Let us now turn to the problem of determining the 
electromagnetic force acting on an atomic system that 
is initially prepared in an arbitrary internal (electronic) 
quantum state. For convenience, we shall employ the mul- 
tipolar formalism. On recalling Eqs. (pij) and to- 
gether with Eq. dini), we find that Eq. (HH) [with E(fA) 
and B(fA) replaced by E'(fA) and B'(fA), respectively] 
can be rewritten as 



dujV 



1 d^ 
-d X 

uv at 



dE'(r,w) 



e.g., Ref. [101). Note that we are free to choose a con- 
venient operator ordering in Eq. (|106p . because E'(r,w) 
commutes with d. 



A. General case 

In order to calculate the average force as a function 
of time, we first formally integrate the Heisenberg equa- 
tions of motion for the fundamental fields i'^{r,uj,t) to 
obtain the source-quantity representation of E'(r,a;,t). 
The result reads (see Appendix IE| 

E'(r,c^,t) =Ej„,(r,w,i) + E^ourcc(r,w,i), (109) 
where 

E^_(r,c.,t)=E'(r,c.)e— ' (110) 

and 

Eourcc(r:W:i) 

f dt'e-^-(*-*')lmG[r,fA(t'),^]d(i')- (m) 



Substituting Eq. (fT09)) together with Eqs. pTO]) and pTTI) 
into Eq. (|106p . we arrive at 

(F(t)) = (Ff,ec(0) + (rsourco(t)), (112) 

where 



J_ d_ 

iuj dt 



(d(t) x[VxE^_(r,..,t)]) 



r=rA(t) 



H.c. 



(113) 



and 



Here, 



VxE'(r,c.)ll +H.C., (106) (f1_(<)) = T dc^c.^ /Ve-^-(*-*') 

J J r=f A i ^ Jo Jo 



where E'(r,w) is defined according to Eq. ([8]). Decom- 
posing F into an average component (F) (where the ex- 
pectation value (...) is taken with respect to the internal 
atomic motion and the medium-assisted electromagnetic 
field only) and a fluctuating component 



we may write 



AF = F - (F), 



(F) + AF. 



(107) 



(108) 



In the following, we will only consider the average force 
(F) (for a discussion of the force fluctuation (AF^), see. 



X V(d(t)Im G[r, f A(t'), Md(t')) \ + H.c. 

(115) 

is the electric part of the average force associated with the 
source-field part of the medium-assisted electromagnetic 
field, and 

(FTot^M = |^^°"d^c.A^*di'e— 

X {d{t) X {V X ImG[r, f A(t'), ^]}d(i')) j + H.c. 

J r=fA(t) 

(116) 
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is the respective magnetic part. Equations (|112p - ||116p 
are still general in the sense that they apply to both 
driven and nondriven atomic systems and to both weak 
and strong atom-field coupling. 



B. Nondriven atom in the weak-coupling regime 

When the atomic system is not driven, i.e., 

(117) 

then (Ff,.oc(i)) =0. Consequently, the average force, re- 
ferred to as CP force, is determined by the source-field 
part only. 



(F(t)) = (F^ourccW). 



(118) 



Even more specifically, we assume that the density op- 
erator of the initial quantum state of the field and the 
internal (electronic) motion of the atomic system reads 



Q=m){{0'}\^a, 



(119) 



where the density operator of the internal motion of the 
atomic system a can be written as 



(120) 



(^nm = with \n), |m) being the internal atomic 

energy eigenstates). In order to calculate the dipole- 
dipole correlation function appearing in Eqs. (|115p and 
(I116p . we make use of the expansion 



(121) 



and write 



{d{t) d{t')) 

= H (inin<Sldm'n'{Amn{t)Am'n'(t')). (122) 



m,n ra'n' 



In the weak-coupling regime, the Markov approxi- 
mation can be exploited and the correlation functions 
{A-mn{t)Am'n'{t')) cau be calculated by means of the 
quantum regression theorem (see, e.g., Ref. [1^). For 
this purpose, the (intra-atomic) master equation has to 
be solved for arbitray initial conditions, which in gen- 
eral requires knowledge of the specific level structure of 
the atomic system under consideration. Only if the rel- 
evant atomic transition frequencies are well separated 
from each other, one can go a step forward constructing 
a general solution. In this case, the off-diagonal density- 
matrix elements can be regarded as being decoupled from 
each other and from the diagonal elements. We find (see 
Appendix [F| 



(^■Arfjijiit^ Arfj^i ^i' {t )^ ^ nm' (^-A^-fi^i* it )^ 



X e 



{iLD„„(fA)-[r™(fA)+r„(fA)]/2}(t-t') 



(123) 



{t>t\ n). Here, 

(f-A) = COnin + SuJra{rA) " SUnirA) (124) 

are the body-induced position-dependent shifted transi- 
tion frequencies [ta =fA(^)], where 

5t^,^(fA)=5]<5cc>f„(fA), (125) 

k 

jdfe^ImG^i) (fA,fA,cj)d, 



with 



Wmfe(fA) - W 



and 



r™(fA) = ^r:;(fA) 



(126) 
(127) 



are the position-dependent level widths, with 

r!;(fA) = ^ep™fe(fA)]p™fe(fA)]' 

xdfemlmG [fA,fA, Wmfe(fA)]d™fc. (128) 

One should point out that the position-independent (in- 
finite) Lamb-shift terms resulting from G*^°^ (f a, f a, ^) 
[recall Eq. ([SU])] have been thought to be absorbed in 
the transitions frequencies w„m- Equation (jl26[) can be 
rewritten by changing to imaginary frequencies [cf. the 
discussion below Eq. (j63p ]. resulting in 

Suj^^ivA) = -■^0[a)„fc(rA)][tI'mfe(rA)]^ 



dfo„ReG^^^ [fA,fA,'imfe(fA)]d„ife 
du u'^ojkmirA) 



/""j..„2~ ^dkTnG^'^\rA,rA,iu)dmk 



TTft Jq 



(129) 



Recall that in the perturbative treatment the vdW po- 
tential of an atomic system in a state |m) is identified 
with the energy shift hSujm, so it is not surprising that 
Eq. together with Eq. corresponds to Eq. §^ 

together with Eqs. §^ and §7]), if in Eq. the <I>„fe 

are replaced with ujmk- The calculation of 



(^mn(^)) — Cnm(i) 



(130) 



[o'nm(O) = cr„m] then leads (under the assumptions made) 
to 



anmit) = e{*'^""(^A)-[r™(fA)+r„(i-A)]/2}t^ 



(131) 



for niy^n [cf. Eq. (|123p ]. so the remaining task consists 
in solving the balance equations 

(Tmm(i) = -r„(fA)fT™„(t) + ^ (f A )fT„„ (t) . (132) 



With these preparations at hand, the CP force can 
be calculated in the following steps. We first substi- 
tute Eq. together with Eqs. p^ and P^ into 
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Eqs. (|115p and (|116p and perform the time derivative in 
Eq. (jll6p . Introducing slowly varying density-matrix ele- 
ments (Tnm{t) = e*'^"'"*(T„m(t), we then perform the time 
integrals in the spirit of the Markov approximation, by 
making the replacements crnm{t') ^ <^nm{t) as well as 
fA(^') '-^ ^A{t) and letting the upper limit of integra- 
tion tend to infinity. Recalling Eq. (I118|) together with 
Eq. (|114p . we derive 



(F(t)) =^(T„™(i)F„„(fA), 



(133) 



This result requires two comments. First, in Eqs. (|135p 
and (|136p the replacement G{r,rp^,uj) i-^ G(^^(r, ta,^) 
has again been made, which can be justified by simi- 
lar arguments as in Sec. IIIII [cf. the discussion preceding 
Eq. Second, from the derivation of Eqs. (fn3| - p36l) 

it is clear that these equations arc valid provided that 
the center-of-mass motion can be regarded as being suf- 
ficiently slow. More precisely, they hold if the condition 

G[r, rA{t+At),Lu] w G[r, taW, cj] for At < T'^ (137) 



F™„(fA) = F^J.„(fA) + F^f (fA), (134) 



where 



V0d™fcImG(i)(r,rA,t^)dfc„ 1 , ^ 

X r 7Z — ^ TTTZ — 7-. — ^ — — — rvTTT > + H.C. 



+ (ife„(fA) - i[rfe(fA) + r„(rA)]/2 J 
and 

d,„fc X [V X ImG(^)(r,fA,^)]dfc„ 
U) + UknirA) - i[rfc(fA) + rm(fA)]/2 



(135) 



H.C. 



is satisfied, where Tq is a characteristic intra-atomic de- 
cay rate. Under this condition, the internal (electronic) 
and external (center-of-mass) motion of the atomic sys- 
tem decouple in the spirit of a Born-Oppenheimer ap- 
proximation. As a result, ta effectively enters the equa- 
tions as a parameter, so that the caret will be removed 
in the following (rA ^ ta). 

Wc finally rewrite Eqs. (|135p and (|136p . by using con- 
tour integration and going over to imaginary frequencies 
[cf. the discussion below Eq. ([63]) ]. RecaUing the defini- 
tion of a.mni'^) = CKmri(rA, 'J-') sls givcu in Eq. (|69p and 
introducing the abbreviating notation 

) = <^nfc(rA) + i[Tm{rA) + rfe(rA)]/2, (138) 

we derive 

Ft,(rA) = F^^^-(rA) + F«^';;(rA), (139) 

Frf(rA) = Fr„^'-(rA)+F— (rA), (140) 



(136) where 



Fmn'(rA) = -J ^ /" duu^[{a„,nh{^A,iu) + {a„,nh{^A,-iu)]VGlf{v,VA,iu)\ , (141) 



Fmn(rA)= <^AioE®('^"fe)"""'t(rA)V(»d„feGW[r,rA,r!™„fe(rA)]dfe„i + H.c. (142) 



and 



27r Jo 



du u^Tt 



Wmn(rA) T / ■ \ ^mn{^A) -J , . n 

■■ "mn(i"A, iw) ■ Q:m„(rA, -m) 



[VxG(i)(r,rA,m)] 



r=rA 

(143) 



F™?f''(rA) = < /ioE®('^"'=)'^™"('"A)f^mnfe(rA)d™fe X (V X G(i)[r,rA,annfe(rA)]dfe„) \ +H.C. (144) 



[(TrT)j ^Tm\. Equation together with Eq. pM)) and Eqs. P^ ^ pii)) is the natural generalization of 
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Eq. §^ together with Eqs. ([65]), §7^1, and The 
above result is the first nonperturbative expression for 
the CP force that incorporates its time dependence in 
case of excited atoms and correctly accounts for body- 
induced shifting and broadening of atomic transition 
lines. 

In the short-time limit, Tc't^ 1, Eq. (|133p reads 

{F{t)) ~ (F(0)) - ^c7,„„(0)F,„„(rA), (145) 

which for (Jnm{^) — ^ni^mi rcduccs to 

(F(t))^(F(0))=F^/(rA). (146) 

For the nonrelativistic Hamiltonian (|46p . we can always 
choose real dipole matrix elements {dmn =d„m), reveal- 
ing that dmn dnm is a symmetric tensor so that, recall- 
ing Eq. (|17p . we may exploit the rule 

S,jVG[f{r,r,uj) = 25,,\iGg^ (s, r, w)|,=r, (147) 

which is valid for any symmetric tensor 5'. Hence, 
Eqs. (fT4T|) and (fT42ll [together with Eq. ^] lead to 

F^;'°^(rA) = r duu'[{ai),,irA,iu) 

47r Jq 

+ {ai)ij{rA,-iu)]VAG'lf{rA,rA,iu) (148) 

and 

Ff/ivA) = ^J2^Mnl{rA){v 

k 

®dikG^^\v,vMvA)]dki] +H.C. (149) 

i r=rA 

[Sl/fc(rA) = il/ifc(rA)]- Ignoring the position-dependent 
shifts and broadenings of the atomic energy levels, i.e., 
disregarding the position dependence of the atomic po- 
larizability [a.i{YA,iu)^af\iu)], Eqs. (fT48| and (fT49)) 
reduce to the perturbative result in Eq. ([68|) together 
with Eqs. (ESI), UnZl), and ^ [F^/(rA) Fi(rA)]. Note 
that this result can be obtained without choosing real 
dipole matrix elements [oniiu) + a.i[—iu) being symmet- 
ric in this case]. In the long-time limit, Fct^ 1, Eq. (|133p 
obviously reduces to ground-state force 

(F(t)) ~ ^a„„(^)F„„(rA) = F'^^ti^f^) (150) 

m,n 

[Foo'(rA) = 0], because of cJnm{co) = 5no5mo- 

As already mentioned, the expression for the ground- 
state CP force Foo(rA) obtained in lowest-order pertur- 
bation theory, Eq. (|77|) . agrees with the expression ob- 
tained from LRT. However, its naive extrapolation in 

the sense of the replacement a.^^\uj) ^ olq{vat^) in 
Eq. ([77)) [25! is wrong, because it results in Eq. (|148p with 
2(ao)ij(rA, iu) instead of (ao)y (rA, w)-|-(ao)ij(rA, -iu). 
As a result, a noticeable influence of the level broadening 



on the off-resonant part of the CP force is erroneously 
predicted in Ref. [25|] (cf. Sec. IV Cp . thus demonstrating 
that body-induced level broadening is a nonperturbative 
effect which lies beyond the scope of the LRT approach 
to the problem. 

Equation (|148p reveals that even the ground-state CP 
force cannot be derived from a potential in the usual way, 
because of the position dependence of the atomic polariz- 
ability. Nevertheless, it is a potential force, provided that 
it is an irrotational vector, i.e., 

Va X Foo(rA) = j Auu^ ^ | [VAWfco(rA)] 

+ [VATk{vA)\-g^\[{ao)^J{vAM 

-(- (Q:o)jj(rA,-m)] X VAG'-j^(rA,rA,m) = 0. (151) 

While for effectively one-dimensional problems (e.g., for 
an atom in the presence of planarly, spherically, or cylin- 
drically multilayered media) this condition is satisfled, 
there are of course situations where it is violated, imply- 
ing that Eq. (|148p is inaccessible to perturbative methods 
in principle. 

When the atomic system is initially prepared in a co- 
herent superposition of states such that cr„m(0) 7^ is 
valid for certain values n and m with n ^ then — 
according to Eq. (|133p — the corresponding off-diagonal 
force components (Jnm(t)Fmn{'^A) can also contribute to 
the total force acting on the atomic system. Interestingly, 
such transient off-diagonal force components contain con- 
tributions not only from the electric part of the Lorentz 
force but also from the magnetic part, as can be eas- 
ily seen from inspection of Eqs. (|143p and (|144p . Thus 
an atomic qubit = (|0) + |1))/V2 (cf., e.g., Ref. [47]) 
near a body feels, in electric dipole approximation, both 
an electric and a magnetic force in general. 

Let us briefly comment on atomic systems displaying 
(quasi)degeneracies, i.e., systems exhibiting transitions 
with bJmn — ^yn'yi' (jn^m! and/or n^n'). In such a case, 
the assumption that the (relevant) off-diagonal density- 
matrix elements decouple from each other as well as from 
the diagonal ones can no longer be made. Let us assume 
that the degenerate sublevels are not connected via elec- 
tric dipole transitions [dmm' =0 if Wmm' —0). The degen- 
eracy related to the different possible projections of the 
angular momentum of an atom (in free space) onto a cho- 
sen direction is a typical example. Taking into account 
that the degeneracy is removed when the atom is close 
to a body, it may be advantageous to change the basis 
within each degenerate sublevel accordingly and consider 
the master equation in the new basis. An equation of the 
form of Eq. (|123p is then valid in the new basis. Note that 
the new basis will in general depend on the position of 
the atom, thus introducing an additional position depen- 
dence of the CP force. While Eq. p3ip also remains valid 
in the new basis for Umn 7^ 0, this is not in general true 
for the temporal evolution of the density-matrix elements 
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with Umm' — SO that, instead of the balance equations 
(|132p . a system of equations has to be solved in which di- 
agonal density-matrix elements and off-diagonal elements 
with ujmm' — are coupled to each other. 



C. Example: Excited atom near an interface 

To illustrate the effects of body-induced level shift- 
ing and broadening, let us consider a two-level atom 
with (real) transition dipole matrix element cIa = dio — 
dp^{cos4>sin9 ex + sin<l)s\n9 Gy + cosO Gz) (doo^dn =0), 
which is situated at position za very close above (z > 0) a 
semi-infinite half space [z < 0) containing a homogeneous 
dispersing and absorbing magnetodielectric medium. Let 
Suj = Suji — Sluq denote the (position-dependent) shift of 
the transition frequency. Using the Green tensor in the 
short-distance limit, from Eqs. and we 

derive (see Appendix [G| 



Suj{za) = (5w,.(za) -I- SlooAza), 



(5w,.(za) 



C |e[wio(zA) 



ftz^ |e[wio(zA) 



(152) 



(153) 



2Cojio{za) 
hnz"^ 



du 



e{iu) — 1 



Cj^^izA) + u"^ e{iu) + 1 ' 

(154) 



where 



C 



did 



327ren 



(155) 



Note that in the short-distance limit the medium effec- 
tively acts like a dielectric one. Since the relation Qiq = 
u!io + Slu is valid, Eq. (|152p together with Eqs. (|153p and 
(|154p is a highly transcendental equation for the deter- 
mination of 6u!. To solve it, we first note that the off- 
resonant term SlOq^ may be neglected in most practical 
situations. For example, for a single-resonance medium 
of Drude-Lorentz type. 
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FIG. 1: (a) Transition frequency shift (solid and dotted lines) 
and (b) decay rate (solid and dotted lines) versus bare transi- 
tion frequency for a two-level atom that is situated at distance 
ZA from a semi-infinite half space medium of complex permit- 
tivity according to Eq. (I156|l and whose transition dipole mo- 
ment is perpendicular to the interface [cijp/tJT = 0.75, '^/lut 
= 0.01; a;|di/(37rfeoc^) = 10"^ za/At = 0.0075 (solid and 
dashed lines), 2a/At =0.009 (dotted and dot-dashed lines)]. 
For comparison, the approximate results obtained by using 
the bare frequencies in Eqs. (|158|) and (|159|l are also displayed 
(dashed and dot-dashed lines). 



e{Lo) 



1 



(156) 



and the parameter values in Fig. [H one can easily verify 
the inequality 



6uJor{zA) 

i^w{za) ~ 2ftz^cj^(I)io(zA) 



< 



< 10" 



(157) 



Thus, keeping only the resonant part of the frequency 
shift, we may set 



Suj{za) 



C \e[[dw{zAW-l 
■ftzi|e[(I-io(zA)]-Hl|2 



(158) 



For e(cl'io) from Eq. (|156p . Eq. (|158p is a fifth-order poly- 



nomial conditional equation for 6uj, which may be solved 
numerically. Having calculated Suj, we may calculate the 
(position-dependent) decay rate r = ri. Neglecting the 
small free-space decay rate, we replace the Green tensor 
by its scattering part as given by Eq. (jG4p . hence from 
Eqs. P?7P and (f^g)) we obtain 



r(^A) = 



4C Im£[cDio(zA)] 
hzl |e[(I'io(zA)] + lp 



(159) 



The resonant part of the CP force on the excited atom 
in the short-distance limit can be found by taking the 
derivative of the scattering part of the Green tensor 
[Eq. (|G4p ] with respect to za and substituting the re- 
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suit into Eq. (fT49l) (Z = 1). We derive (F'^^ = Fl^^e^) 
3C |£[Oio(za)]P-1 

^11 (^A = r uro ( M I 112 ' (160) 

zX \e[ilio[ZA)\ + Ir 
where, according to Eq. (|138p . 

r!io(zA) = tDio(zA) + jr(zA)/2. (161) 
Using Eq. p56p . we see that (7,r<Cti'T) 

^2 

e[fiio(2A)] = 1 + ^ Tj-T — N — 7- — rr— 7 — r ■ 

- '^io(^a) - 4r(^A) + tJ'^ioI^^a) 

(162) 

Equation (|160p differs from the perturbative result in 
two respects. First, the bare atomic transition frequency 
LUiQ is replaced with the (position-dependent) shifted fre- 
quency uJiQ. Second, the absorption parameter 7 of the 
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FIG. 2: The resonant part of the CP force Fi^X^ x 10-V(3C) 
on a two-level atom that is situated at distance (a) za/Xt ~ 
0.0075 and (b) za/At ~ 0.009 of a semi-infinite half space 
medium of complex permittivity according to Eq. (|156p and 
whose transition dipole moment is perpendicular to the inter- 
face (solid lines). The parameters are the same as in Fig. [T] 
For comparison, both the perturbative result (dashed lines) 
and the separate effects of level shifting (dotted lines) and 
level broadening (dash-dotted lines) are shown. 



medium is replaced with the sum of 7 and the (position- 
dependent) atomic decay rate F. The sum 7 -I- F obviously 
plays the role of the total absorption parameter. 

The dependence of Suj and F on loiq in the short- 
distance limit is shown in Figs. [TJa) and (b), respec- 
tively, and Fig. [2] displays the resonant part of the CP 
force as a function of luiq. From Fig. [2] it is seen that in 
the vicinity of the (surface-plasmon induced) frequency 
ujs = \/ uj^ + ujp/2 an enhanced force is observed, which 
is attractive (repulsive) for red (blue) detuned atomic 
transition frequencies ujiq < (wio > ws) — a result al- 
ready known from perturbation theory (dashed curves in 
the figure). However, it is also seen that due to body- 
induced level shifting and broadening the absolute value 
of the force can be noticeably reduced (solid curves in 
the figure). Interestingly, the positions of the extrema of 
the force remain nearly unchanged, because level shifting 
and broadening give rise to competing effects that almost 
cancel. 

In order to calculate the off-resonant part of the CP 
force on the excited atom in the short-distance limit, we 
first note that, according to Eq. 



'^io(za) 



4dA' 



'Cjfo{zA) + [u + T{zA)/2]^ 

, Cjf,{zA)+u^ + r^{zA)/4 
•c^2^(za) + [u-F(za)/2]2 



(163) 



Substituing Eq. (|163|) into Eq. (|148p and making use of 
Eq. (|G7|) [where f{u) is given by times Eq. p63p ]. we 
derive (Ffi =F{'[e^) 



^^iT(^a) 



3C 



du 



enu) 



1 



e{iu) + 1 
wio(^a) 



' CoI^{za) + [u + T{za)/2Y 
^u:Uza)+u^+T\za)/^ 



Colo{zA) + [u-T{zA)/2f 
Note that for a two-level atom the relation 



FS^{za) = -FiT(za) 



(164) 



(165) 



is valid. 

In Fig. [31 the off- resonant part of the CP force is shown 
as a function of the bare atomic transition frequency. 
Obviously, the shift of the transition frequency has the 
effect of raising and lowering the perturbative values of 
the force (dashed curves) for wio <a;s and loiq >ujs, re- 
spectively, which is in full agreement with the frequency 
response of the frequency shift shown in Fig. [IJa). The 
influence of the decay rate on the CP force is extremely 
weak, as it can be seen from the insets in the figure. 
This may be understood by the fact that in contrast to 
the case of the resonant part of the CP force, where the 
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FIG. 3: The off-resonant part of the CP force F^iX^ x 10"^ 
/(3C) on a two- level atom that is situated at distance (a) 
za/Xt = 0.0075 and (b) za/Xt = 0.009 of a semi-infinite half 
space medium of complex permittivity according to Eq. (|156|) 
and whose transition dipole moment is perpendicular to the 
interface (solid lines). The parameters are the same as in 
Fig. [T] For comparison, the perturbative result (dashed lines) 
is shown. The insets display the difference between the force 
with and without consideration of the level broadening (solid 
lines). For comparison, we show this difference in the case 
where the level shifts are ignored (dashed lines). 



decay rate enters directly via the Green tensor, the in- 
fluence on the ofF-resonant part is more indirect via the 
atomic polarizabihty. Due to the specific dependence on 
the atomic polarizabihty, the leading-order dependence 
is quadratic in T and not linear in T as erroneously pre- 
dicted from LRT [lEl- Physically, the weak influence of 
the level broadening on the off-resonant part of the CP 
force may be regarded as being a consequence of the fact 
that this part corresponds to energy nonconserving pro- 
cesses (the energy denominators being nonzero), which 
implies that they happen on (extremely short) time scales 
where real photon emission does not play a role. 

Comparing the magnitudes of the resonant and off- 
resonant components of the CP force, we see that the 
off-resonant component is smaller than the resonant one 



by about two orders of magnitude. However, this obser- 
vation should be considered with great care. While the 
two-level atom is a good model for calculating the reso- 
nant part of an atom in an excited state, such a simplifi- 
cation is not justified in general when all higher levels can 
contribute to the off-resonant force component. However, 
provided that the convergence of the corresponding sum 
is sufficiently fast, we can still conclude that the resonant 
part of the CP force is dominant. 



VI. SUMMARY 

Basing on electromagnetic-field quantization that al- 
lows for the presence of dispersing and absorbing lin- 
ear media, and starting with the Lorentz force acting 
on a neutral atom, we have extended the concept of CP 
force beyond the well-known results derived on the ba- 
sis of normal-mode quantization or LRT in leading order 
of pertubation theory to allow for (i) magnetodielectric 
bodies, (ii) an atom that is initially prepared in an arbi- 
trary internal (electronic) quantum state, thereby being 
subjected to a time-dependent force, (iii) the position de- 
pendence of the force via the atomic response, and (iv) 
arbitrary strength of the atom-field coupling. The basic 
formulas also apply to the calculation of the radiation 
forces arising from excited fields such as the force acting 
on a driven atom. 

For a first analysis, we have restricted our attention to 
a nondriven atom in the weak-coupling regime, so that 
the internal atomic dynamics can be treated in Markov 
approximation. It turns out that the force is a superposi- 
tion of force components weighted by the time-dependent 
intra-atomic density-matrix elements that solve the intra- 
atomic master equation. Each force component is ex- 
pressed in terms of the Green tensor of the electromag- 
netic field and the atomic polarizabihty, which — through 
the position-dependent energy level shifts and broaden- 
ings — now depends on the position of the atomic system. 
In consequence even the force components resulting from 
the electric part of the Lorentz force cannot be derived 
from potentials in the usual way. Clearly, the position de- 
pendence of the atomic polarizabihty become noticeable 
only for very small atom-body separations. In order to 
illustrate the effect, we have considered a two-level atom 
in the vicinity of a planar semi-infinite medium. 

When the atomic system is initially prepared in an 
eigenstate of its internal Hamiltonian, then only force 
components associated with diagonal density-matrix ele- 
ments appear. They solely result from the electric part of 
the Lorentz force and reduce to the CP forces obtained 
in lowest-order perturbation theory if the atomic polariz- 
abihty is replaced with its position-independent pertur- 
bative expression. Force components that are associated 
with excited intra-atomic energy levels are of course tran- 
sient. As in the course of time an initially excited level 
is depopulated and lower lying levels are populated, the 
force that initially acts on the atomic system in the ex- 
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cited state changes with time to the force that acts on 
the atomic system in the ground state. 

The results further show that when the atomic system 
is initially prepared in an intra-atomic quantum state 
that is a coherent superposition of energy eigenstates, 
then additional force components associated with the cor- 
responding off-diagonal density-matrix elements are ob- 
served. Thus an atomic qubit would typically feel such 
off-diagonal force components. It should be pointed out 
that not only the electric but also the magnetic part 
of the Lorentz force can contribute to the off-diagonal 
force components, with the magnetic contributions be- 
ing proportional to the transition frequencies. Clearly, 
off-diagonal force components are transient. 

In contrast to the transient force components that are 
associated with excited energy levels, off-diagonal force 
components carry an additional harmonic time depen- 
dence. Clearly, if the oscillations are too fast, it can be 
difficult to detect them experimentally, since they may 
effectively average to zero. In this case it may be advis- 
able to assign them to the fluctuating part of the force 
rather than to the average force. The situation may be 
different in cases where strong atom-field coupling (not 
considered here) gives rise to Rabi oscillations. 
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APPENDIX A: DERIVATION OF THE 
MULTIPOLAR HAMILTONIAN (l38l) 



To perform transformations of the type 

O' = Jjdu\ (Al) 

with tj being given by Eq. together with Eq. ([37]), 
we apply the operator identity 



1 



2! 



(A2) 



Recalling the commutation relations (HJ and ([5]) , it is not 
difficult to prove that the basic fields f (r, lo) are trans- 
formed as 

f;(r,c.)=fA(r,c.) + -i- /'dVPi(r')G*(r',r,^). (A3) 



Using Eq. (|A2p together with the commutation relation 



[eo£^fe(r),A;(r)]=:iM^;(r-r'), cf. Ref. [33 , we find that 



E'(r)=E(r) + lpi(r). 



(A4) 



To transform the momenta of the charged particles, the 
identities 



Va(5(r- f-A - Ar^) 



A-1 



\5. 



ruA 



a/3 



VJ(r-fA- AF/3), (A5) 



dAf„V(5(r-fA-AFc«) = '5(r-fA) - <5(r-f„), (A6) 



dAArQ,V5(r — f a — Ar^) 



-6{r - 



Ca) + / dA^(r - f-A - Af-a) 
Jo 



(A7) 



are helpful. They can be proved with the aid of the def- 
initions (j28p and (|29|). and via (partial) integration with 
respect to A. Using Eqs. (|A5|) - (|A7p we derive 



Pa = Pa - <ZaA(f„) - / dVH„(r) X B(r), (A8) 



where 3a (r) is defined as in Eq. ([39]) . Further, the fol- 
lowing quantities remain unchanged under the transfor- 
mation (jAip . because they commute with both A(r) (cf. 
Ref. ^7}) and fa, 

A'(r) = A(r), B'(r) = B(r), ^'(r) = ^(r), (A9) 

f^=r„, f-A^r-A, PA(i") = PA(r), ^aW = 'y5A(r), (AlO) 

P:,(r)^PA(r), 0'Jr)-0a(r), HUr) = H„(r). (All) 

Applying the transformation rules (|A3|) . (|A8p . and 
(IA9p - (|Alip . we may now express the minimal-coupling 
Hamiltonian ([1]) in terms of the transformed variables. 
Recahing Eq. ([HI) together with Eqs. ([8|)-([10|) and mak- 
ing use of the relations and 



/ da;4lmG(r,r',c.)-5<5(r-r') 
(cf. Ref. [3^), we derive 



(A12) 



/dVP:,^(r)P:,^(r) 



d'^rP^(r)E'^(r) 
n 2 



d'rp'Air)ip'Air)+ / dV p:,(r)^'(r). (A13) 



In order to simplify the last two terms of Eq. (|A13p , we 
recaU Eq. §^ as well as P^(r) = eoV(^^(r) and E'll(r) 
= — Vi^'(r), obtaining with the aid of partial integration 



i / d3rP:,(r)V^:,(r)+ / dV P:,(r) V^'(r) 



^ d\P'Ar)P'l{v)- I d\P'A{r)I^'Hr). (A14) 
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Combining Eqs. (|A13p and (jA14|l . and noting that inte- 
grals containing mixed products of transverse and lon- 
gitudinal vector fields vanish, we obtain Eq. ((38|) . where 
we have made use of Eqs. (|A10I) and (jAlip and hence 
dropped the primes of all quantities containing the par- 
ticle coordinates only. 

In the simpler case in which the center-of-mass coor- 
dinate is treated as a parameter, the transformation law 
[8|) changes to 

Pa = Pa - 9aA(f,) ~ / d^Q^^ X B(r). (A15) 



Equations (|A3|, dAH, and (lA9 |) - (|Alip remain formally 
the same, provided that the replacement f ai— is made. 



APPENDIX B: ORDERS OF MAGNITUDE OF 
INTERACTION TERMS 

To estimate the order of magnitude of atom-field in- 
teractions, let us introduce the typical atomic length and 
energy scales 



Ob 



^0 = Z^gE^ 



Z2^13.6eV 



(Bl) 
(B2) 



(aB, Bohr radius; £'r, Rydberg energy), where me and 
— e are the electron mass and charge, respectively, Z^se 
is the typical effective nucleus charge felt by the electrons 
giving the main contributions to the interaction terms to 
be calculated, and ao = e^/(47reo^c) is the fine-structure 
constant. As a rough estimate we can then make the re- 
placements 



qa e, Too 
— > ao, ta — 



TOe, ujki Eo/n, (B3) 
p = meEoao/h (B4) 



-(') 

V, Pa 



[for the last replacement, see Eq. (jC7p ]. With regard to 
the length scale of variation of the medium-assisted elec- 
tromagnetic field we may make the replacements 

V ^ X^^^uj/c, Vip Vip^ujA, (B5) 
eW ^ Er^ujA, BW ^ B-(cj/c)A (B6) 

(A*^'^ A). Noting that materials typically become 
transparent for frequencies that are greater than 20 eV 
(cf. Ref. HI), 

£(r, lu)kI^ G^i) (r, r', tj) w for huj > 20 eV, (B7) 

we should require that 



hw < 20 eV 



< 1. 



(B8) 



With these approximations at hand, the orders of mag- 
nitude of Ai_E defined by Eq. ([55]) and A2E defined by 
Eq. ((54|) in Sec. IIII Al can be estimated to be 

„2/l2 p2 2/12 

AliJ ^ 1^ = ^-^Eo = g'Eo = O(g') (B9) 



and 



2to„ 



1 



-Bo + ^ 



eV^^ eaoS/LpepA 2 2^ 2 
5 h 2 1- e aQVf 



1 + 2 



E'n 



fku 
E'o 



En 



1 + huj/Eo 



where the dimensionless coupling constant 
g = eaoA/h 



(BIO) 
(Bll) 



has been introduced. Note that in Eq. (jBlOp we have 
approximated Pa P, because in Sec. IIIII we treat an 
atom at rest, hence relative and absolute momenta are 
identical. 

In order to give a rough idea of the magnitude of the 
coupling constant g, we need to estimate the magnitude 
of the field strength A. In the context of the current work 
we consider interactions of an atomic system with the 
vacuum electromagnetic field, so the relevant quantity is 
the vacuum fluctuation of the field strength. Recalling 
Eqs. ([8]) and pTjl and making use of the commutation 
relations ([4]) and (O as well as the integral relation p9|) . 
we find 

([AE(rA)]^) = ({0}|E2(rA)|{0}) - ({0}|E(rA)|{0})2 

fe roo 12 

= / dc^'— ImTrG(rA,rA,c^')- (B12) 

TTEO Jo C 

When the atomic system is placed sufficiently far away 
from all macroscopic bodies, a good estimate for the in- 
tegral can be given by using the vacuum Green tensor 
ImG^*'^ (r, r, w) = ti;/(67rc)J, leading to 



([AE(rA)] = 



(B13) 



where cj is a characteristic frequency contributing to the 
interaction, cf. Eq. (jB8p . Hence making the replacement 



A^ 

[cf. Eq. dm], we find 
g ^ ZcS\ 



"0 



10' 



(B14) 



(B15) 



depending on the specific atomic system considered and 
the characteristic frequencies of the medium. When the 
atom is situated close to some macroscopic body, the 
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scattering Green tensor becomes much larger than the 
vacuum Green tensor, and the approximation leading to 
Eq. I|B13|1 is not valid anymore. The increased value of 
the coupling constant g is reflected by the failure of the 
perturbative result for small atom-surface separations. 

The orders of magnitude of the contributions of the 
three terms in Eq. ([5^ to the eigenvalue shift in Sec. lIII Bl 
can be estimated according to 



|dE'(rA)r 



(eaoE) 



Eq 



1 + huj/Eo 
-O(g^) 



1 + huj/Eo 



(B16) 



2m 



h{ujki + uj) 

2 



En 



eaopB \ . 
2me / 1 + hu/Eo 

Eq 

l + hw/Eo 



0[{Zcsao9f], 



(B17) 



2 1 



Eo = 0[{Z,fiaogy]. (B18) 



Next, let us estimate the orders of magnitude of 
the various contributions to the Lorentz force given in 
Sec. lIV Al The magnitudes of the first and third terms in 
curly brackets in Eq. (|97p can be approximated according 
to 



|dE(fA)| ~ eaoE = g{huj)^0{g), 



(B19) 



fABW(fA) X d 

- ^eaoBv = i (^5) {hu) = 0{gv/c). (B20) 



In order to estimate the magnitude of the second term, 
we make use of the relation 

maYa = Pq - gQA(fa) (B21) 

in order to introduce relative momenta [recall Eq. 
leading to 



^^f„B(rA) X 



E 



2ma 



PaB(fA) X ia 



+ ^F«B(fA) X A(f a) + ifAB(f a) X d. 



2ma 



(B22) 



Combining this with 



7r^PaB(fA) X Va 

= (Zoffaog) = 0{Z^saog), (B23) 



2me 



E^r„B(fA)xA(fA; 



2nia 

e^aoAB 
2mp 



and Eq. (|B20|) . we see that the magnitude of the sec- 
ond term in curly brackets in Eq. (|97p is 0{ZaSCiog + 
Zcsott^g^ + gv/c) — 0[{Zcsao + v/c)g]. The magnitudes 
of the different contributions to Eq. (|98|) are 



d X B(fA) 

epB e 



E ^ [P" ~ <?aA(f a)] X B(fA) 



^AB 



= .9(1 + 25) 



= 0(g), 



(B25) 



d X B(fA) 



d X 



eaoLuB = 5 ( ^) ) = 0(g), (B26) 



JaVa (g> B(fA) -f B(f a) ® VaJ-a 



eaQVLoB 



(- 

Vc 



(B27) 



Finally, let us compare the contributions of the Ront- 
gen interaction to the temporal evolution t\(r,Lu,t) with 
that from the electric dipole interaction. 



^t^{t)d{t) X (Va X G*[rA(<),r,w]; 



j-dit)Gl[iAit),v,co] 



(B28) 



X B[rA{t),t]d{t) X (Va x G^[rA(i), r, u;]) 



0[iZ,eaofg]. 



rUe 

(B29) 
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APPENDIX C: CALCULATION OF THE 
PERTURBATIVE CORRECTIONS ^ AND (USD 

Recalling Eq. ^ together with Eqs. (|8))- (fT0)) . making 
use of the commutation relations ^ and (O, and apply- 
ing Eq. HI]), Eq. ^ leads to 



a A—e.7n 

where we have introduced the notation 

^(ll)G-^(IIHr,r',o.) 

= j ^""-^ j d3s'<5^(lini--s)G(s,s',cj)J^(ll)(s'-r'). 

(C2) 

Applying the sum rule 

y^I^T^ry^^uidiu^du+dki^diu), (C3) 

OL k 

we can rewrite Eq. (jCip as 



/"OO 

AiEi = —y / dwwfeidz/cIm-^G-^(rA,rA,a;)dfc;. 
k ^0 

(C4) 

To calculate A2-E, as given by Eq. (|54p . we first calculate 
the matrix elements therein. Recalling Eqs. (|4|)- ([T0|) . we 
obtain 



(/|({0}|dV^(r)r=rJ{lA(r,c^)})|fc) 

= -d;fctlGA(rA,r,c^), (C5) 

m„ 



-ai({0}|^^p„A(rA)|{lA(r,c.)})|fc) 



dik G\{rA,r,u), 

UJ 



(C6) 



where the second matrix element has been obtained by 
means of the identity 



{l\Pa\k) = -iukidik- 



(C7) 



Substituting Eqs. ((C5)) and (|C6|l into Eq. (JM]), we then 
derive 



J„ ujki+Lu J 



d^ridlk)^ 



k A— e,m 



+ ^(^G.),„^(rA,r,c.)(^G^)^„(rA,r,^y 
-E^ r^^d,4-^'lmllGll(rA,rA,c.) 

+ ojkiLo [im "G^(rA, TA, t^) + Im ^G'I (fa, ta, t^) 
-cc>2;Im^G^(rA,rA,c^)ldfc,, (C8) 



where we have again made use of the identity (jl9|l . 
Adding Eqs. ((C4| and ((C8)) according to Eq. ([ItI), on 
using the identity G^^G^ + ^G" + I^G^ + "GH [which 
directly follows from the definition l|C2p together with 
S{r) =5ll(r) + (S-*- (r)] , we eventually arrive at Eq. ([55]) . 

The derivation of Eq. (|75|) is completely analogous. The 
relevant matrix elements can be calculated with the aid 
of Eq. (HI]) together with Eqs. (©-(dU]) and the commuta- 
tion relations (|4]) and (O, cf. the remarks below Eq. ([40]). 
The result is 

-(/|({0'}|dE'(rA)|{lUr,^)})|fc) = -d;feGA(rA,r,c.). 

(C9) 

Substituting Eq. ([09]) into Eq. ([75]) yields 

dcij 



d^r 



k X—e.m 

X {dik)i{dki)3 {Gx)^^{yk, r, ^) {Gl) (fa, r, w), 



(CIO) 



from which Eq. ([75]) follows by means of Eq. (fTO)) . 



APPENDIX D: EQUIVALENC E OF LORENTZ 
FORCES dMl) AND (fT05]) 



To transform the first term in Eq. (|105p . we apply the 
the rule (jA4p . recall that integrals over mixed products 
of transverse and longitudinal vector fields vanish, and 
use the identity for the first term in Eq. (|A14[) as well as 
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Eqs. dm) and (fT03)) . We thus derive 



the commutation relations ([4]) and we obtain 



VA^dV [PA(r)E'(r)] 

= VA^dV [PA(r)E(r)] + ^VAydV [PA(r)Pi(r)] 
= Va^ dV [PA(r)E(r)] + ^-Va ^ dVPi(r) 

-Va^ d^rpA(r)(^A(r) 
= Va /dVPA(r)E(r). (Dl) 



In order to simplify the second term in Eq. ()105|) . we use 
the definitions IMl), ([Ml), and (Ei) to calculate 



Q 

= ^^1a[@a{r) X fa - f-Q X (r) 

a 

~ ^ X! ^" X f A - f A X 0a(r) 

+ i[PA(r) X i-A -i-A X PaW 
= MA(r)+MR(r). 



(D2) 



Consequently, recalling that B'(r) =B(r), we may write 

Va / dV i ^ [H„(r) X f„ - f„ X H„(r)] B'(r) 
= Va y d\ [MA(r) + MR(r)] B(r) (D3) 



as well as 



d_ 

di 



dVPA(r) X B'(r) 



d 
dt 



dVPA(r) X B(r) 



(D4) 

Substituting Eqs. (|Dl|), jDSl, and (|D4} into Eq. pUS)) . 
we see that Eq. (|105p is equivalent to Eq. 



APPENDIX E: EQUATIONS OF MOTION FOR 



In electric dipole approximation, the temporal evolu- 
tion of the basic fields f^(r, uj, t) is governed by the Hamil- 
tonian given in Eq. together with Eqs. (|45p . (HS)) . and 
(jSn)) . Using Eqs. ([5]) and (I1T])-(1221) (with the unprimed 
fields being replaced with the primed ones) and applying 



= -tujt'xir, + \ dWGA[f aW, r, lo] 



2huji 



fAWd(t) X (Va X G^[fA(t),r,H 



^—\dit)xB'[rj,{t),t]dit) 



x(VAxG^[rA(t),r,c.])}. (El) 



The third and fourth terms in Eq. (|E1|) , which are due to 
the Rontgen interaction, are smaller than the second one 
by factors of v/c and ^(Zoffao)^, respectively [Eqs. (|B28P 
and (|B29[) in Appendix |B] , so according to the nonrela- 
tivistic approximation, Eq. (jEl[) reduces to 



f;(r, t) = ^iuji'^{v, ^.t) + - d{t)Gl[vA{t), r, (E2) 

which can be integrated to yield [f^(r, io, 0) = f^(r, ljj)] 

fA(r,a;,i) = fj,f,ee(r, tj, t) + f;,o^,„^(r, w, t), (E3) 
where 

f;ft,,,(r,c.,i)-e— *f;(r,c.) (E4) 

and 

f;„(r,c.,t) = -^y'di'e— (*-*')d(^')G^[fA(^'),r,H• 
(E5) 

Substituting Eqs. (1E3|-(1E5| into Eq. © [E(r,u;,i) i-^ 
E'(r,Li;,i)] and using the identity p9|) . we arrive at 
Eqs. (frogi) - (|TTTD . 



APPENDIX F: INTRA-ATOMIC EQUATIONS OF 
MOTION 



An estimation similar to that given for the fields 
f'(r,tt',t) shows that in the nonrelativistic limit the sec- 
ond term in the interaction Hamiltonian in electric dipole 
approximation (|50p can be disregarded in the calculation 
of the temporal evolution of the intra-atomic operators 
Amn{t). By representing the (unperturbed) intra-atomic 
Hamiltonian in the form of Eq. (|53|) . recalling Eqs. ([8]) 
and (|12ip . and applying standard commutation relations, 
it is not difficult to prove that the A„m (t) obey the equa- 
tions of motion 



Amn — ~r\H , An 



A 



Akn) / dcjE'(fA,t^) 



-/ dcjE'^(fA,t^)(d„fei„fe-dfe,„ifc„) 
Jo 



(Fl) 
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We now substitute the source-quantity representation for 
E'(rA,w) — E'[rA(i), w, (and its Hermitian conjugate) 
according to Eqs. (fT09l) - (|TTT|) into Eq. (|Fl]). Carrying 
out the time integral in the source-field part in Eq. (|F1[) 
in the Markov approximation, we may set, on regarding 
ta — ^A{t) as being slowly varying. 



dt^E^ourco(l"A,w) ^^g,nn{rA)A„^n, (F2) 



where 



.(fA) 



/ dcj oj^Im G(rA, f A) t^)d,i 

Jo 



X C['^nm(fA)-'^] 



(F3) 



[C{x) = t:5{x) + iV/x], with ujnmi^A) being the shif- 
ted transition frequencies. Substituting Eq. (|F2|) into 
Eq. jFT]), we obtain 



Ar. 

with 



+B. 



h ^ 

k 



i 

h 



y^.'i-nkgmlAlk 



k,l 



X! dfc„g^;ii„ (F5) 



[rn^n)^ and 



fefc 



(F6) 



with 



= ^ |dmfcgfeiAm/ — dfcmg^;^;,„| 

X! (dfcmgm/^fei - d„fcg;„;i/fc|, (F7) 
k,l^k 

where F^n denotes contributions from the free-field part 
in Eq. (|109p . Taking expectation values with respect 
to the internal atomic motion and the medium- assisted 
electromagnetic field, with the density-matrix given by 
Eq. ()119p . we can use the property (|117p . finding that 
the terms -Fmn do not contribute. In the absence of (qua- 
si)degeneracies such that 



(F8) 



we may disregard couplings between different off-dia- 
gonal transitions and between off-diagonal and diagonal 
transitions and thus omit the terms Bmn and Cmm , hence 
upon using the decomposition 



(F9) 



where Suj^{rA) and rJ;(fA), respectively, are defined ac- 
cording to Eqs. (fT26| and (fT28| [with G^^^ (f a, f a, t^) in- 
stead of G(fA,fA,t^) in Eq. fm]) ]. Eqs. ^ and (|F6] 
lead to Eqs. (pg)) . (fT?T|) . and p^ . 



APPENDIX G: HALF SPACE MEDIUM 

The equal-position scattering Green tensor for a semi- 
infinite half space which contains for z < a homo- 
geneous, dispersing, and absorbing magnetodielectric 
medium reads for z > Esl 



(F4) G(i)(r,r,c.) = 



dq J- e'^''"^ 
Po 




where 



-/32 
-/32 
2^2 



ePo-f3 
£/3o + /3 



(Gl) 



(G2) 



are the reflection coefficients for s- and p-polarized waves, 
respectively {i3q = uj'^ / — q"^ with lm/?o>0, (3'^ = e fj.uj'^ / 
~ q^ with Im/3 > 0). For q > \ijj\/c and q > -^/jipjlcjl/c, 
respectively, the approximations 



/3o ^ «g, /? =i iq 



(G3) 



can be made. Due to the exponential factor the integra- 
tion interval is effectively limited to values q<l/ z. In the 
short-distance limit z-\/|£/7[|a;|/c<C 1, we therefore intro- 
duce a small error, if we extrapolate the approximations 
?3|) to the whole integral, resulting in 



e{u:) - 1 




(G4) 



Note that the magnetic properties of the medium rep- 
resented by the permeability /i begin to contribute via 
terms proportional to 1/z. Substitution of Eq. (jG4p 
into the first term of Eq. p29p for Sujnk = Sujio yields 
Eq. (fT55D . 

In order to obtain Eq. (I154p . we recall Eq. (|Gip to 
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write 



dw/(w)G(i)(r,r,m) 

/'OO nOO 

^ / du/(u) / d6oi 



1 
X < I 1 





P 9 



6g 
bl 



2&2-(H)2 



(G5) 



having changed the integration variable to the imaginary 
part of /3o {Po=ibo). Let ujm be a characteristic frequency 
of the medium such that 



£(iu) — 1 <C 1 for w > ujm- 



(G6) 



For u > lom, the approximation /3 ~ /3o holds, and conse- 
quently the reflection coefficients Tj, rp are independent 



of bo. The frequency integral effectively extends up to 
frequencies of the order c/z, hence in the short-range 
limit zlum/c <^ 1 c/z 3> wm) we introduce only a 
small error by extrapolating this approximation to the 
whole frequency integral. Performing the bo integral, re- 
taining only leading-order terms in uz/c (in consistency 
with zlum/c ^1) we derive 



duf{u)G'^^\r,r,iu) 



327rz3 



J{u) ejiu) - 1 
e{iu) + 1 




Using Eq. jGTl [with f{u) = 
Eq. (fT29l) . we obtain Eq. p54)) . 



(G7) 
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